O 

(N 

o 

Q 

(N 




Axial Nucleon form factors from lattice QCD 



P. 

(a) 



C. Alexandrou (a < b) , M. Brinet W, J. Carbonell W, M. Constantinou (a \ 
A. Harraud W, P. Guichon K. Jansen e , T. Korzec (aJ) , M. Papinutto (c) 

Department of Physics, University of Cyprus, P.O. Box 20537, 1678 Nicosia, Cyprus 
O Computation-based Science and Technology Research Center, 
Cyprus Institute, 20 Kavafi Str., Nicosia 2121, Cyprus 
' c ' Laboratoire de Physique Subatomique et Cosmologie, 
UJF/CNRS/IN2P3, 53 avenue des Martyrs, 38026 Grenoble, France 
^ CEA-Saclay, 1 RFU '/Service de Physique Nucleaire, 91191 Gif-sur-Yvette, France 
(e) NIC, DESY, Platanenallee 6, D-15738 Zeuthen, Germany 
^' Institut fur Physik Humboldt Universitat zu Berlin, Neuitonstrasse 15, 12489 Berlin, Germany 



(N 
> 

in 

oo 

o 

(N 



We present results on the nucleon axial form factors within lattice QCD using two flavors of 
degenerate twisted mass fermions. Volume effects are examined using simulations at two volumes 
of spatial length L — 2.1 fm and L = 2.8 fm. Cut-off effects are investigated using three different 
values of the lattice spacings, namely a = 0.089 fm, a — 0.070 fm and a = 0.056 fm. The nucleon 
axial charge is obtained in the continuum limit and chirally extrapolated to the physical pion mass 
enabling comparison with experiment. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 



I. INTRODUCTION 



The nucleon (N) form factors are fundamental 
hadronic observables that probe the structure of the nu- 
cleon. Experiments to measure the electromagnetic nu- 
cleon form factors have been carried out since the 50's. 
A new generation of experiments using polarized beams 
and targets are currently under way at major facilities 
in order to measure the nucleon form factors more accu- 
rately and at higher values of the momentum transfer. 
The nucleon form factors connected to the axial vector 
current are more difficult to measure and therefore less 
accurately known than its electromagnetic form factors. 
As in the electromagnetic case the nucleon matrix ele- 
ment of the axial vector current is written in terms of 
two Lorenz invariant form factors, the axial form fac- 
tor Ga{<1 2 ) and the induced pseudo-scalar form factor, 
G p (q 2 ) where q 2 is the momentum transfer squared. The 
nucleon axial charge gA — Ga(0), which can be deter- 
mined from /3— decay, is known to a high precision. The 
q 2 -dependence of G^q 2 ) has been studied from neutrino 
scattering [l[ and pion electroproduction 2, 3] processes. 
The nucleon induced pseudo-scalar form factor, G p (q 2 ), 
is even less well known. Muon capture at low q 2 val- 
ues [H and pion electroproduction for larger Q 2 @, Q are 
the main experimental sources of information for G P (Q 2 ). 
Both G^g 2 ) and G p (q 2 ) have been discussed within chi- 
ral effective theories 0, Q ■ I n this work we present results 
on these form factors obtained in lattice QCD using two 
degenerate light quarks (Np=2) in the twisted mass for- 



mulation Q- 

Twisted mass fermions Q provide an attractive for- 
mulation of lattice QCD that allows for automatic 0(a) 
improvement, infrared regularization of small eigenval- 
ues and fast dynamical simulations. For the calculation 
of the nucleon form factors in which we are interested 
in this work, the automatic 0{a) improvement is particu- 
larly relevant since it is achieved by tuning only one pa- 
rameter in the action, requiring no further improvements 
on the operator level. 

The action for two degenerate flavors of quarks in 
twisted mass QCD is given by 

S = S g + a 4 ^ x(x) [D w +m clit +i~f 5 T 3 n] x(x) , (1) 

X 

where D\y is the Wilson Dirac operator and we use the 
tree-level Symanzik improved gauge action S g [9[. The 
quark fields \ are m the so-called "twisted basis" ob- 
tained from the "physical basis" at maximal twist by a 
simple transformation: 
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[1 + ir 3 -f 5 ]x and vb- 



1 r q i 



(2) 



We note that, in the continuum, this action is equivalent 
to the standard QCD action. A crucial advantage is the 
fact that by tuning a single parameter, namely the bare 
untwisted quark mass to its critical value m cr , a wide 
class of physical observables are automatically 0(a) im- 
proved [7| . A disadvantage is the explicit flavor symme- 
try breaking. In a recent paper we have checked that this 
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breaking is small for the baryon observables under consid- 
eration in this work and for the lattice spacings that we 
use To extract the nucleon FFs we need to eval- 

uate the nucleon matrix elements (N(p', s')\A^\N(p, s)), 
where \N(p', s')), \N(p,s)) are nucleon states with fi- 
nal (initial) momentum p'(p) and spin s'(s). Due to its 
isovector nature, the axial vector current, defined by 



(3) 



receives contributions only from the connected diagram 
for a = 1, 2 and up to 0(a) for a = 3. Simulations includ- 
ing a dynamical strange quark are also available within 
the twisted mass formulation. Comparison of the nu- 
cleon mass obtained with two dynamical flavors and the 
nucleon mass including a dynamical strange quark has 
shown negligible dependence on the dynamical strange 
quark [l5|. We therefore expect the results on the nu- 
cleon form factors to show little sensitivity on a dynam- 
ical strange quark as well. 

The axial current matrix element of the nucleon 
(N(p r , s')\A"(0)\N(p, s)) can be expressed in terms of the 
form factors Ga and G p as 
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(N(p',s')\Al\N(p,s))=i 



N 



E N (p')E N (p) 



u N (p',s') 



g^75 

2tojv 



G P (q 2 



-u N (p,s). (4) 



In this work we consider simulations at three values 
of the coupling constant spanning lattice spacings from 
about 0.05 fm to 0.09 fm. This enables us to obtain re- 
sults in the continuum limit. We find that cut-off effects 
are small for this range of lattice spacings. We also ex- 
amine finite size effects on the axial form factors by com- 
paring results on two lattices of spatial length L = 2.1 fm 
and L = 2.8 fm [IM3- 



II. LATTICE EVALUATION 



A. Correlation functions 

The nucleon interpolating field in the physical basis is 
given by 

J(x) = e abc [u aT (x)C l5 d b (x)] u c {x) (5) 

and can be written in the twisted basis at maximal twist 
as 

J(x)=^=[t + ij 5 ]e abc \u aT (x)Cj 5 d b {x)\ u c {x). (6) 

The transformation of the axial vector current, A^(x), 
to the twisted basis leaves the form of A^ (x) unchanged. 
The axial renormalization constant Za is determined 
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FIG. 1: Connected nucleon three-point function. 



non-perturbatively in the RI'-MOM scheme using two 
approaches [ltl [2(| and [U [22[ both of which yield con- 
sistent values. We use the values of Z A found in latter 
approach [2l[ , which employs a momentum source l23ll 
and a perturbative subtraction of 0(a 2 ) terms (241. |25||. 
This subtracts the leading cut-off effects yielding only 
a very weak dependence of Z A on (ap) 2 for which the 
(ap) 2 — > limit can be reliably taken. It was also shown 
with high accuracy that the quark mass dependence of 
Za is negligible. We find the values 



Z A =0- 757(3), 0.776(3), 0.789(3) 



(7) 



at (3—3.9, 4.05 and 4.2 respectively. These are the values 
of Za which we use in this work to renormalize the lattice 
matrix element. 

In order to increase the overlap with the proton state 
and decrease overlap with excited states we use Gaussian 
smeared quark fields [26J, |27| for the construction of the 
interpolating fields: 

y 

F = (l + aH) n , 

3 

H(x,y:U(t)) = ^[ U ii x ) S ",v-t + U K X ~ *)<Wd • 

i=l 

In addition, we apply APE-smearing to the gauge fields 
Up entering the hopping matrix H. The smearing pa- 
rameters are the same as those used for our calculation 
of baryon masses with a and n optimized for the nucleon 
ground state [l3j]. The values are: a = 4.0 and n = 50, 
70 and 90 for (3 = 3.9, 4.05 and 4.2 respectively. 

In order to calculate the nucleon matrix element of 
Eq. ([J]) we calculate the two-point and three-point func- 
tions defined by 

G{q, t f ) e-**'-*rg Q (J a (t f ,S f )J fi {U,Si)) (9) 



where T° and T k are the projection matrices: 

r° = i(i + 7o ) ) T k = iT° l5lk . 



(10) 
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FIG. 2: The ratio of Eq. (|11[> for representative momentum 
combinations at /3 = 3.9 and different values of fi. The 
filled (black) circles show results with a sink-source separa- 
tion tf/a = 14 and the filled (red) squares for tf/a = 12, 
shifted to the left by one time-slice. 



The kinematical setup that we used is illustrated in 
Fig. Q] We create the nucleon at ti=0, at Xi=0 (source) 
and annihilate it at a later time t f with ft =0 (sink) . The 
current couples to a quark at an intermediate time t. 
For our kinematics q = —ft. From now on all quanti- 
ties are given in Euclidean space and where Q 2 = —q 2 is 
the Euclidean momentum transfer squared. The leading 
time dependence of the Euclidean time evolution and the 
overlap factors cancel for the ratio 



R»(T,q,t) 



G" (I\ q, t) G{fti, t f - t)G(0, t)G(0, t f ) 



(11) 

yielding a time-independent value 



lim lim 

tf — t— >oo t — 1»— K 



(12) 



We refer to the range of t- values where this asymptotic 
behavior is observed within our statistical precision as 
the plateau range. As mentioned already, only the con- 
nected diagram contributes. It is calculated by perform- 
ing sequential inversions through the sink yielding the 



form factors at all possible momentum transfers and cur- 
rent orientation /x. Since we use sequential inversions 
through the sink we need to fix the sink-source separa- 
tion. Statistical errors increase rapidly as we increase 
the sink-source separation. Therefore we need to choose 
the smallest possible that still ensures that the nucleon 
ground state dominates when measurements are made at 
different values of t. In order to check that a sink-source 
time separation of ~ 1 fm is sufficient for the isolation 
of the nucleon ground state we compare the results at 
(3=3.9 obtained with tf/a=12 i.e. tt ~ 1 fm with those 
obtained when we increase to tf/a=14 [ll|. As can be 
seen in Fig. [21 where we plot i? M (T, q, t), the plateau val- 
ues are compatible yielding the same plateau value for 
the two time separations. This means that the shorter 
sink-source separation is sufficient and the ground state 
of the nucleon dominates in the plateau region. We there- 
fore use in all of our analysis tf — t{ ~ lfm. 

New inversions are necessary every time a different 
choice of the projection matrix T a is made. In this work, 
we consider choices, which are optimal for the form fac- 
tors considered. Namely we use the spatial T's and con- 
sider the spatial component of the current i.e. we extract 
the form factors from 



TT(T\q)-- 



4m 



qkqi 

2m 



G p (Q 2 )-(E+m)S hk G A {Q 2 



(13) 



where k = 1,2, 3, and c = 



2 m 2 



E(E+m) ■ 



B. Simulation details 

The input parameters of the calculation, namely /3, 
L/a and a/x are summarized in Table HI The lattice spac- 
ing a is taken from the nucleon mass as described in the 
next section. The pion mass values, spanning a mass 
range from 260 to 470 MeV, are taken from Ref. [Hj]. At 
TOtt 300 MeV and /3=3.9 we have simulations for lat- 
tices of spatial size L=2.1 fm and L=2.8 fm allowing to 
investigate finite size effects. Finite lattice spacing effects 
are studied using three sets of results at /3=3.9, /3=4.05 
and /3=4.2 for the lowest and largest pion mass avail- 
able in this work. These sets of gauge ensembles allow us 
to estimate all the systematic errors in order to produce 
reliable predictions for the nucleon axial form factors. 



C. Determination of the lattice spacing 

The nucleon mass has been computed on the same en- 
sembles that are now used here for the computation of the 



nucleon axial form factors 13( . Therefore we can use the 
nucleon mass at the physical point to set the scale. We 
show in Fig. |31 results at three values of the lattice spac- 
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P = 3.9, a = 0.089(1)(5) fm, r /a = 5.22(2) 



24 X 48, L = 2.1 rm 
32 3 x 64, L = 2.8 fm 


_ r\ r\r\Ar\ r\ r\c\£i a r\ nnoc a riin 

a/x 0.0040 0.0064 0.0085 0.010 
Stat. 943 553 365 477 

(GeV) 0.3032(16) 0.3770(9) 0.4319(12) 0.4675(12) 
m n L 3.27 4.06 4.66 5.04 
a\i 0.003 0.004 
Stat. 667 351 
m n (GeV) 0.2600(9) 0.2978(6) 
rrivrL 3.74 4.28 


,3 = 4.05, a = 0.070(1)(4) fm, r /a = 6.61(3) 


32 3 x 64, L = 2.13 fm 


a^i 0.0030 0.0060 0.0080 
Stat. 447 326 419 
m n (GeV) 0.2925(18) 0.4035(18) 0.4653(15) 
m*L 3.32 4.58 5.28 


/3 = 4.2, o = 0.056(1)(4) fm r /a = 8.31 


32 3 x 64, L = 2.39 fm 
48 3 x 96, L = 2.39 fm 


ay, 0.0065 
Stat. 357 

(GeV) 0.4698(18) 
m,L 4.24 
ay 0.002 
Stat. 245 
m-n (GeV) 0.2622(11) 
m n L 3.55 



TABLE I: Input parameters (P, L,ay) of our lattice calculation and corresponding lattice spacing (a) and pion mass (ra^). 



ings corresponding to /3=3.9, /3=4.05 and /3=4.2. As can 
be seen, cut-off effects are negligible and we can therefore 
use continuum chiral perturbation theory to extrapolate 
to the physical point. For the observables discussed in 
this work the nucleon mass at the physical point is the 
most appropriate quantity to set the scale. This also 
provides a cross-check for the determination of the lat- 
tice spacing as compared to the pion decay constant used 
in the meson sector. If lattice artifacts are under control 
then these two determinations should be consistent, un- 
der the assumption that quenching effects due to the ab- 
sence of the strange quark from the sea are small for these 
quantities. In order to correct for volume effects we use 
chiral perturbation theory to take into account volume 
corrections coming from pions propagating around the 
lattice, following Ref. [29[ . A similar analysis was carried 
out in Ref. [ill at /3 = 3.9 and (3 — 4.05 and we refer to 
this publication for additional details. In addition, this 
work is extended by an analysis of results at j3 = 4.2. In 
Table [TTI we give the volume corrected nucleon mass. 

To chirally extrapolate we use the well-established 
C(p 3 ) result of heavy baryon chiral perturbation theory 
(HBxPT) given by 

3<7 2 

m N =m%- 4 Cl ml - -f^m 3 . (14) 

We perform a fit to the volume corrected results at /3=3.9, 
^=4.05 and ^=4.2 and extract r =0.462(5) fm. Fit- 



am n 




aniN 


amjv(L — >• oo) 


P = 3.9 


0.2100(5) 


5.04 


0.5973(43) 


0.5952 


0.1940(5) 


4.66 


0.5786(67) 


0.5760 


0.1684(2) 


4.06 


0.5514(49) 


0.5468 


0.1362(7) 


3.27 


0.5111(58) 


0.5043 


0.1338(2) 


4.28 


0.5126(46) 


0.5115 


0.1168(3) 


3.74 


0.4958(43) 


0.4944 






P = 4.05 




0.1651(5) 


5.28 


0.4714(31) 


0.4702 


0.1432(6) 


4.58 


0.4444(47) 


0.4426 


0.1038(6) 


3.32 


0.4091(60) 


0.4056 






P = 4.2 




0.1326(5) 


4.24 


0.380(3) 


0.3763 


0.0740(3) 


3.55 


0.306(4) 


0.3049 



TABLE II: Results on the nucleon mass. The last column 
gives the values after a volume correction. 



ting instead to the /3=3.9 and /3=4.05 results we find 
ro=0.465(6) fm showing that indeed cut-off effects are 
small. To estimate the error due to the chiral ex- 
trapolation we use HBxPT to C(p 4 ), which leads to 
ro=0.489(ll). We take the difference between the 0(p 3 ) 
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FIG. 3: Nucleon mass in units of ro at three lattice spacings 
and spatial lattice size L such that m^L > 3.3. The solid 
(black) and dashed (red) lines are fits to 0(p 3 ) and 0(p 4 ) 
HBxPT. The physical point is shown with the asterisks. Re- 
sults at P = 3.9 and 24 3 x 48 are shown with filled (red) 
circles, at /3 = 3.9 and 32 3 x 64 with the filled (blue) squares, 
at /3 = 4.05 and 32 3 x 64 with the filled (green) triangles, at 
P = 4.2 and 32 3 x 64 with the open (yellow) square and at 
P = 4.2 and 48 3 x 96 with the star (magenta). 



and 0(p 4 ) mean values as an estimate of the uncertainty 
due to the chiral extrapolation. Fits to other higher 
order x?T formulae are also shown in Fig. [3] These 
are described in Ref. [l3[ and are consistent with 0(p 4 ) 
HBxPT. Using r =0.462(5)(27) and the computed r /a 
ratios we obtain 

a^3. 9 = 0.089(1)(5), 
0/3=4.05 = 0.070(1)(4), 
a^=4.2 = 0.056(2)(3). 

These values are consistent with the lattice spacings de- 
termined from /„■ and will be used for converting to phys- 
ical units in what follows. We note that results on the nu- 
cleon mass using twisted mass fermions agree with those 
obtained using other lattice 0(a 2 ) formulations for lat- 
tice spacings below 0.1 fm [H| . 



III. RESULTS 

In the first subsection we discuss our results on the nu- 
cleon axial charge and in the second subsection we discuss 
the momentum dependence of the axial Ga(Q 2 ) and the 
induced pseudo-scalar G P (Q 2 ). 



A. Axial charge 

Our lattice results on the nucleon axial charge are 
shown in Fig. [4] and listed in Table IIIII In the same 
figure we also show results obtained using Np = 2 + 1 




FIG. 4: The nucleon axial charge. Results using Nf = 2 
twisted mass fermions are shown using the same notation as 
that of Fig. [3] Crosses show results obtained using Nf = 2 + 1 
DWF, circles are results in a mixed action approach on a 
lattice of size 20 3 x 64 and the triangle on a lattice of size 
28 3 x 64. 



domain wall fermions (DWF) by the RBC-UKQCD col- 
laborations [3(| and using a mixed action with 2+1 fla- 
vors of asqtad sea and domain wall valence fermions by 
LHPC 31]. The first observation is that results at our 
three different lattice spacings are within error bars. The 
second observation is that results at two different volumes 
are also consistent. The third observation is that there 
is agreement among lattice results using different lattice 
actions even before taking the continuum and infinite vol- 
ume limit. 



rriTr 


Lm-K 


9A 


gA(L -> oo) 


P = 3.9 


0.4675 


5.04 


1.163(18) 


1.167 


0.4319 


4.66 


1.134(25) 


1.140 


0.3770 


4.06 


1.140(27) 


1.150 


0.3032 


3.27 


1.111(34) 


1.133 


0.2978 


4.28 


1.103(32) 


1.106 


0.2600 


3.74 


1.156(47) 


1.162 


P = 4.05 


0.4653 


5.28 


1.173(24) 


1.177 


0.4035 


4.58 


1.175(31) 


1.182 


0.2925 


3.32 


1.194(66) 


1.218 


P = 4.2 


0.4698 


4.24 


1.130(26) 


1.144 


0.2622 


3.55 


1.138(43) 


1.146 



TABLE III: Results using Nf = 2 twisted mass fermions 
(TMF) on the axial nucleon charge. The last column gives 
the values after a volume correction. 
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FIG. 5: The nucleon axial charge as a function of Lm, 
notation is the same as that of Fig. U 
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Finite volume effects 



In order to assess volume effects we plot in Fig. [5] re- 
sults on gA versus Lm^. Besides TMF results we show 
the results obtained using Np = 2+1 DWF [3(j as well as 
within the mixed action approach |3lj |. As can be seen 
the results are consistent with a constant in the whole 
range of Lm, spanned. In particular we do not observe 
a decrease in the value of gA for values of Lm, ~ 3.3. 
Therefore, given that finite volume effects are negligible 
for the smallest value of m,L = 3.3 as compared to the 
value we find at m^L = 4.3, we conclude that for all of 
our data for which m^L > 3.3 volume effects are small. 



We can estimate the volume correction on gA within 
heavy baryon chiral perturbation theory (HByPT) in the 
so called small scale expansion (SSE) 32], which includes 
explicitly the A degree of freedom. In this scheme one 
expands the results in powers of a small parameter e, 
which denotes small pion four-momenta, the pion mass, 
baryon three-momenta and the nuclcon-A mass splitting 
in the chiral limit. Writing 



gA(m n , oo) = 3a(tott, L) - Ag A (m^,L) (15) 
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FIG. 6: The nucleon axial charge as a function of the lattice 
spacing in units of ro at: rom-n = 0.615, roro T = 0.85, rom n = 
0.95, rom T = 1.10, from top to bottom. We use ro = 0.462(5) 
extracted from the nucleon mass. In the upper most and 
lower most graphs we show both the linear (dotted line) and 
constant fits (dashed line). 



is given by |33j 
(9a) 3 K 



9a 



°ml ^ Ki(L\ri\m*) 



E 



6^/2 
4 f2f> 



E 



Ki{L\n\m n ) 



K (L\n\m w ) ..... 

i 

dyy^2 K o(L\n\f(m W) y)) 

ft 

L\H\f(m^,y) 

Ai(L|n|/(m ff ,y)) 



2 _o r°° 



277T2/2 



d ^E^— 



K (L\n\f(m w ,y)) 



Ki{L\n\f{m^,y)) 



the dependence of gA on the spatial length L of the lattice 



'27tt/2 A ^ L\n\u, 



L\n\f(mir,y) 
y — ) e- i|fi|m - + 0(e 4 ) (16) 
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with /(tojt, y) = yjm\ + y 2 + 2yA and f n the pion de- 
cay constant in the chiral limit which we approximate 
with its physical value i.e. we take / w = 0.092 GeV. In 
the sum Y^h an vectors n are summed except n = 0. In 
order to estimate the volume correction AgA we take the 
experimental value of the axial charge in the chiral limit 
i.e. g A ~ g A xp — 1.267 and the nucleon - A mass splitting 
in the chiral limit Ao = 0.2711. For the A axial coupling 
constant we use the SU(4) relation g% = 9g A * p /5 and for 
the nucleon to A axial coupling constants ca — 1.5. The 
estimated volume corrected gA is given in Table IIII1 

In order to assess cut-off effects we use the simula- 
tions at three lattice spacings at the smallest and largest 
pion mass used in this work. We take as reference 



pion mass the one computed on the finest lattice and 
interpolate results at the other two /3-values to these 
two reference masses. In Fig. [5] we show the value of 
gA at these reference pion masses computed in units of 
tq. We perform a fit to these data using a linear form 
g A {a) = <?a(0) + c(a/ro) 2 . The resulting fit is shown in 
Fig. HI Setting c = 0we obtained the constant line also 
shown in the figure. As can be seen, for both large and 
small pion masses the slope is consistent with zero yield- 
ing a value in the continuum limit in agreement with the 
constant fit. Therefore we conclude that finite a effects 
are negligible and for the intermediate pion masses we 
obtained the values in the continuum by fitting our data 
at /3 = 3.9 and (3 = 4.05 to a constant. 





9a{P = 3.9) 


g A (/3 = 4.05) 


9a{P = 4.2) 


g A (a -)■ 0) 


g A (L — > oo, a — > 0) 


1.1019 


1.165(18) 


1.173(25) 


1.130(26) 


1.159(13) [1.127(40)] 
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1.175(31) 




1.148(21) 


1.155(21) 


0.85 


1.138(28) 


1.179(37) 




1.153(22) 


1.165(22) 


0.686 


1.110(39) 


1.194(66) 




1.127(34) 


1.129(34) 


0.615 


1.153(47) 


1.199(69) 


1.138(43) 


1.154(29) [1.142(76)] 


1.165(29) [1.156(76)] 



TABLE IV: In the second, third and fourth column we give the interpolated values of gA at the value of m^ro given in the first 
column. We used ro/a = 5.22(2), 6.61(3) and 8.31(5) for f3 — 3.9, 4.05 and 4.2 respectively. In the fifth column we give the 
value of gA after extrapolating to a — using a constant fit. In the parenthesis we give the corresponding values when using a 
linear fit. In the last column we give the continuum value of gA for the volume-corrected data. 



The values for g A found at six reference pion masses 
are given in Table [IVJ We give both the continuum values 
obtained using a constant fit when no volume corrections 
are carried out as well with the volume corrected data. 
The volume corrected data extrapolated to a — are 
plotted in Fig. 



2. Chiral extrapolation 



Our simulations cover a range of pion masses from 
about 470 MeV down to about 260 MeV. The pion mass 
dependence for the nucleon axial charge has been stud- 
ied within HBxPT in the SSE formulation [H]. We use 
the one-loop result including explicitly the A degrees of 
freedom in order to extrapolate our lattice results to the 
physical point. We make a three parameter fit to the 
form 
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(18) 



The three parameters to fit are g° A , the value of the 
axial charge at the chiral point, the A axial coupling 
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FIG. 7: The nucleon axial charge obtained by taking the con- 
tinuum limit of the volume corrected data. The shaded area 
shows the best fit to the data shown on the graph. The dashed 
line shows the best fit to the raw lattice data at the three val- 
ues of /3 with the dotted lines being the associate error band. 



constant g\ and a counter-term CsseW- We again take 
the nucleon to A axial coupling constant ca = 1-5, the 
mass splitting between the A and the nucleon at the 
chiral limit, A = 0.2711 and A = 1 GeV. Fitting the 
volume corrected continuum results we find 3^=1-12(8), 
,gi=2.37(1.52) and C SSE = - 1.01(2.01). The parame- 
ters g\ and Csse are highly correlated explaining the 
resulting large error band. Fitting the lattice without 
any volume correction and without extrapolating to the 
continuum limit we obtained £u=l-08(8), gi=2.02(1.21) 
and Csse= — 0.63(1.57) which are consistent with the 
continuum volume corrected results. This shows that 
both cut-off and volume artifacts are small as compared 
to the uncertainty due to the chiral extrapolation. 



B. Axial form factors 

In this section we discuss the results obtained for the 
axial form factors Ga(Q 2 ) and G P (Q 2 ). 

To assess cut-off effects we compare in Fig.[5]results for 
Ga{Q 2 ) and Gp(Q 2 ) versus Q 2 for three different lattice 
spacings at a similar pion mass of about 470 MeV. As 
can be seen, results at these three lattice spacings are 
consistent indicating that cut-off effects are negligible for 
these lattice spacings. We perform a dipole fit to Ga(Q 2 ) 
using 



G A (Q 2 h 



9A 



(l+Q 2 /m 2 A ) 2 



(19) 



with a momentum upper range of Q 2 ~ 1.5GeV 2 . The 
axial mass tua of the fits is larger than in experimen- 
tal value of m e A xp = 1.1 GeV extracted from the best 
dipole fit to the electroproduction data. This is evident 
from the smaller slope shown by the lattice data both for 
twisted mass fermions and domain wall fermions. Assum- 
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FIG. 8 

at TO.„- 



The nucleon axial form factors Ga{Q 2 ) and G P (Q 2 ) 
- 470 MeV at = 3.9 (filled red circles), 4.05 (filled 



green triangles) and 4.2 (yellow squares) versus Q . The line 



is the result of a dipole fit (to the form given in Eq. (|201 
Ga(Q 2 ) (G P (Q 2 )) data on the coarse lattice. 



ing the partially conserved axial current relation and pion 
pole dominance we can relate the form factor G P (Q 2 ) to 
G A (Q 2 ): 



G p {Q 2 )=Ga{Q 2 \ 



G P (0) 
2 2 + m 2 



(20) 



The dependence of these form factors on the pion mass 
is seen in Fig. [3J where we show both Ga and G p com- 
puted at several values of the pion mass spanning a range 
from about 470 MeV to 300 MeV at /3 = 3.9. We show fits 
to the lattice data using a dipole form as given in Eq. (fT9"|) 
for G A {Q 2 ) and to the form given in Eq. (gD| for G P (Q 2 ), 
which described the data rather well. Although the mass 
dependence is weak and the general trend is to approach 
experiment, lattice data show a weaker Q 2 dependence 
as compared to experiment. As already pointed out, the 
best dipole fit to the electroproduction data yields an 
axial mass m e A p = 1.1 GeV [a], and it is shown by the 
solid line. The experimental line in the case of Gp(Q 2 ) 
shown in Fig. [S] is obtained using Eq. (j2U)) and pion mass 
of m p = 135 MeV. 

In Fig. [10] we check for finite volume effects by com- 
paring results obtained at j3 = 3.9 on a lattice of spatial 
length L = 2.8 fm and L = 2.1 fm at m n ~ 300 MeV. 
As can be seen volume effects are negligible for Ga- In 
the case of G p we have a strong dependence on Q 2 as 



9 




FIG. 9: The nucleon axial form factors Ga(Q 2 ) and G P (Q 2 ) 
at /3 = 3.9 for m n = 468 MeV (crosses), m w = 432 MeV 
(filled red circles) and m-n- = 303 MeV (filled blue triangles) 
versus Q 2 . The dashed lines are the result of a dipole fit for 
Ga{Q 2 ) and to the form given in Eq. JUJ for G P {Q 2 ) on the 
coarse lattice. 



Q 2 — » because of the pion pole dependence expected 
for this form factor. Therefore the fits are strongly de- 
pendent on the lowest values of Q 2 that are available. 
E.g. discarding the point at the lowest momentum yield 
the dotted lines which are steeper as compared to includ- 
ing it. Although there is an overall consistency between 
the two data sets at j3 = 3.9 deviations are seen in the 
fits when the same momentum range is used. The fit 
using the whole range of data obtained on the smaller 
lattice, shown by the dashed red line exhibit a weaker 
dependence as compared to the fit using the results on 
the larger lattice but discarding the point at the lowest 
momentum, shown by the dotted blue line. In Figs. [TT1 
and [T2] we compare our results with the results by LHPC 
which were obtained in a mixed action approach that 
uses DWF on staggered sea quarks [31| on a lattice with 
L = 3.5 fm. The results are in agreement in the case of 
Ga{Q 2 ), while in the case of G P {Q 2 ) there are larger dis- 
crepancies. Given the mass dependence of G P (Q 2 ) shown 
in Fig. [!|] a difference in the pion mass of 50 MeV cannot 
fully account for this. Such discrepancies may indicate 
that volume effects are not negligible on form factors such 
as G P (Q 2 ) which are strongly affected by the pion-pole. 

In the case of G P (Q 2 ) one can extract the fit parame- 
ters G p (0) and m p by either fitting the ratio of G p /Ga 



FIG. 10: The nucleon axial form factors Ga and G p at m-n ~ 
300 MeV for a lattice of size 24 3 x 48 (filled red circles) and 
32 3 x 64 (filled blue squares). For Ga(Q 2 ) (upper graph) the 
dotted lines are the best dipole fits to the lattice data. For 
G P (Q 2 ) (lower graph) the dotted lines are fits of lattice results 

to the form ^pLj_^i j discarding the point at the lowest value 

of Q 2 . The dashed lines are fits using all data points. 




0.0 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 

Q 2 (GeV 2 ) 

FIG. If: Axial form factor Ga(Q 2 ) as a function of Q 2 . 
N F = 2 TMF results at m n = 298 MeV are shown with filled 
squares, Nf = 2 + 1 DWF with the crosses at m, = 330 Mev 
and Nf = 2 + 1 using a mixed action of DWF and staggered 
sea quarks at = 356 MeV with triangles. The solid line is 
the dipole parametrization of experimental data. 
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FIG. 12: Induced pseudo-scalar form factor G P (Q 2 ) as a func- 
tion of Q 2 . N F = 2 TMF results at m„ = 298 MeV are 
shown with filled squares and Nf = 2 + 1 using a mixed ac- 
tion of DWF and staggered sea quarks at m n = 356 MeV 
are shown with triangles. The solid line is obtained using the 
parametrization of experimental electroproduction results for 
Ga and pion-pole dominance. 



afi 


m*. (GeV) 


9A 


m A (GeV) 


G,(0) 


m p (GeV) 






P 


= 3.9 






0.0100 


0.4675 


1.163(14) 


1.776(48) 


6.99(74) 


0.738(77) 


0.0085 


0.4319 


1.140(11) 


1.634(32) 


4.40(62) 


0.458(106) 


0.0064 


0.3770 


1.081(31) 


2.021(17) 


3.33(40) 


0.254(113) 


0.004 


0.3032 


1.135(34) 


1.572(82) 


4.30(76) 


0.512(133) 


0.004 


0.2978 


1.160(37) 


1.513(70) 


4.28(63) 


0.459(119) 


0.003 


0.2600 


1.166(28) 


1.445(51) 


2.80(25) 


0.255(95) 








= 4.05 






0.008 


0.4653 


1.174(24) 


1.696(24) 


5.88(43) 


0.595(53) 


0.006 


0.4035 


1.174(32) 


1.725(65) 


4.90(84) 


0.514(124) 


0.003 


0.2925 


1.180(67) 


1.392(44) 


4.13(85) 


0.458(153) 






P 


= 4.2 






0.0065 


0.4698 


1.120(21) 


2.030(93) 


5.72(40) 


0.599(60) 


0.002 


0.2622 


1.158(22) 


1.575(42) 


3.56(18) 


0.325(39) 



TABLE V: Results on the axial nucleon charge and axial mass 
extracted by fitting Ga(Q 2 ) to a dipole form. The two last 
column give the G p (0) and the mass m p by fitting G P (Q 2 ) to 
the form given in Eq. (|20p . 



or use the fitted values for Ga{<1 2 ) and fit to the form 
of Eq. (|20p. The values extracted by performing these 
fits are compatible within error bars. Our lattice data 
on G p (q 2 ) are flatter than pion-pole dominance predicts 
requiring a larger pole mass m p than the pion mass mea- 



sured on the lattice. In Table M we tabulate the resulting 
fitting parameters for all /3 and \i values. The parameters 
Gp(0) and m p have been extracted from fits to the form 
given in Eq. (|20|) . The full set of our lattice results on 
G A {q 2 ) and G P {Q 2 ) is given in the Tables Vl-VIIIof the 
Appendix. 

IV. CONCLUSIONS 

Using Np> = 2 twisted mass fermions we obtain ac- 
curate results on the axial Ga(Q 2 ), and G P (Q 2 ) form 
factors as a function of Q 2 for pion masses in the range 
of about 260-470 MeV. The general feature is a natter de- 
pendence on Q 2 than experiment. This is a feature also 
seen in the electromagnetic nucleon form factors. Finite 
volume effects are found to be small on Ga- Our results 
are in agreement with recent results obtained using other 
lattice fermions such as dynamical Np = 2 + 1 domain 
wall fermions. Having results at three lattice spacings 
enables us to take the continuum limit. We find that 
cut-off effects are small for the values of the lattice spac- 
ings used in this work. Performing a chiral extrapolation 
of our continuum results for the nucleon axial charge, we 
find at the physical point the value qa = 1.12(8). This 
is one standard deviation lower than the physical value. 
The large error associated with our determination of gA 
is mostly due to the chiral extrapolation. Therefore it is 
crucial to perform an analysis with a pion mass closer to 
its physical value. 
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TABLE VI: Results on the axial nucleon form factors at /3 
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TABLE VII: Results on the axial nucleon form factors at /3 = 4.05 
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0.046 


0.759(28) 


3.845(417) 




0.060 


0.734(48) 


3.212(358) 




0.072 


0.634(30) 


2.295(224) 




0.085 


0.584(37) 


1.604(188) 




0.108(1) 


0.440(48) 


1.260(217) 




0.119(1) 


0.414(39) 


0.946(151) 




0.129(1) 


0.364(65) 


0.604(181) 




0.139(1) 


0.328(76) 


0.452(201) 



TABLE VIII: Results on the axial nucleon form factors at j3 



